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Abstract 

Lifespan is influenced by a large number of conserved proteins and gene-regulatory pathways. Here, we introduce a 
strategy for systennatically finding such longevity factors in Saccharomyces cerevisiae and scoring the genetic interactions 
(epistasis) among these factors. Specifically, we developed an automated competition-based assay for chronological 
lifespan, defined as stationary-phase survival of yeast populations, and used it to phenotype over 5,600 single- or double- 
gene l<nockouts at unprecedented quantitative resolution. We found that 14% of the viable yeast mutant strains were 
affected in their stationary-phase survival; the extent of true-positive chronological lifespan factors was estimated by 
accounting for the effects of culture aeration and adaptive regrowth. We show that lifespan extension by dietary restriction 
depends on the Swrl histone-exchange complex and that a functional link between autophagy and the lipid-homeostasis 
factor Arvl has an impact on cellular lifespan. Importantly, we describe the first genetic interaction network based on aging 
phenotypes, which successfully recapitulated the core-autophagy machinery and confirmed a role of the human tumor 
suppressor PTEN homologue in yeast lifespan and phosphatidylinositol phosphate metabolism. Our quantitative analysis of 
longevity factors and their genetic interactions provides insights into the gene-network interactions of aging cells. 



Citation: Garay E, Campos SE, Gonzalez de la Cruz J, Gaspar AP, Jinich A, et al. (2014} High-Resolution Profiling of Stationary-Phase Survival Reveals Yeast 
Longevity Factors and Their Genetic Interactions. PLoS Genet 10(2): e1004168. doi:10.1371/journal.pgen. 1004168 

Editor: Stuart K. Kim, Stanford University Medical Center, United States of America 

Received May 2, 2013; Accepted December 23, 2013; Published February 27, 2014 

Copyright: © 2014 Garay et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This worl( was funded by the Consejo Nacional de Ciencia y Tecnologia (CONACYT CB-1 64889), the University of California Institute for Mexico and the 
United States (UC-MEXUS CN-1 1-559), and a pilot grant from the Northwestern University Physical Sciences-Oncology Center. The funders had no role in study 
design, data collection and analysis, decision to publish, or preparation of the manuscript. 

Competing interests: The authors have declared that no competing interests exist. 

* E-mail: adeluna@langebio.cinvestav.mx 



Introduction 

Lifespan is a complex trait with clear genetic underpinnings 
[1]. Mutations in the insulin/IGF- 1 receptor gene daf-2 increase 
the lifespan of Caenorhabditis ekgans [2] and Drosophila melanogaster 
[3], while mice lacking the insuhn receptor in adipose tissue live 
longer as well [4]. Studies on the budding yeast Saccharomyces 
cerevisiae have also been pivotal to discover conserved longevity 
mechanisms [5-7] . Aging of dividing cells has been modeled in 
yeast by its replicative lifespan, defined as the number of daughter 
cells produced before death [8]. In addition, chronological 
lifespan (CLS), the duration of yeast survival in stationary phase, 
provides a model of aging of non-dividing quiescent cells [9,10]. 
Evidence to date indicates that both forms of yeast aging are 
modulated by the Tor/Sch9 [11,12] and Ras/cAMP/PKA [13] 
pathways. Likewise, the downstream activity of proteins with 
roles in autophagy [14], redox balance [15,16], and mitochon- 
drial function [17,18] influence lifespan from yeast to mammals. 
Moreover, dietary restriction delays both replicative and chro- 
nological aging in yeast; this non-genetic intervention has been 
linked to the down-regulation of Tor and Sch9 in response to 
nutrients [12,19]. 

Current phenotyping strategies allow the characterization of 
yeast CLS in a high-throughput manner. Survival as a function of 



time in stationary phase has been estimated by monitoring the 
outgrowth of either individual-strain cultures in parallel [19,20] or 
large pools of gene-deletion strains coupled to microarray [21,22] 
or deep-sequencing analysis [23]. Such large-scale approaches 
have recently opened the door to the first genome-wide screens for 
lifespan factors [19,21-23]. With these catalogs in hand, the next 
great challenge is to describe how different genetic factors are 
integrated to determine lifespan in yeast. 

Quantitative genetic interaction (GI) analysis has great potential 
to shed light on the genetic architecture of the lifespan phenotype 
and the mechanisms that underlie the aging process. By depicting 
the ways in which individual mutations affect each other's 
phenotypic consequences, GIs inform whether genes act in 
common or independent pathways [24]. A number of break- 
through studies based on the large-scale GI analysis of cell- 
proliferation phenotypes (mainly colony growth) have provided 
important insights into the genetic organization of the living cell 
[25-31]. Multi-conditional and multi-phenotype GI maps have 
further helped to elucidate cellular pathways and responses to 
stimuli [32,33] and, notably, GI networks retain conserved 
features and overall levels of pathway crosstalk across large 
evolutionary distances [34]. Yet, systematic double-mutant anal- 
yses of lifespan phenotypes and their ensuing GI networks are still 
missing. 
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Author Summary 

The budding yeast Saccharomyces cerevisiae has emerged 
as an important model for the genetic analysis of aging, 
and Insights gained about this process in yeast cells 
enhance our understanding of aging in other organisms, 
including humans. Even In yeast, our knowledge of the 
number and Identity of the genes that determine lifespan 
is limited and we are still lacking a general picture of how 
different genetic aging factors work together. Here, using 
an innovative sensitive technique, we have characterized 
the stationary-phase survival of yeast single- and double- 
gene knockout mutants to screen for genes that control 
chronological aging and to score the genetic interactions 
among these genes. Our results showed that an important 
fraction (14%) of the genome contributes to the regulation 
of lifespan, including genes and pathways that had not 
been previously associated to this phenotype. We also 
constructed a genetic interaction map which recapitulated 
llfespan-determlning pathways and highlighted genetic 
associations between the autophagy machinery and the 
phosphatldylinositol-phosphate and lipid-homeostasis 
pathways. Our study provides not only an accurate 
catalogue of yeast aging genes, but also a picture of the 
gene-network wiring of aging cells. 



In this study, we sought to investigate in a quantitative manner 
the genes and gene-network wiring that underUes yeast CLS. By 
using a novel profiling tool based on the competitive survival of 
yeast populations on stationary phase, we generated a high- 
resolution CLS catalogue of 3,878 single-knockout mutations and 
constructed the first systematic lifespan GI map derived from over 
1,800 double-knockout phenotypes. In doing so, we revealed novel 
genetic aging factors and shed light on functional associations 
within and between lifespan-determining pathways. 

Results 

A sensitive high-throughput strategy to profile yeast 
chronological lifespan 

To gain deeper knowledge on which genes shape the lifespan 
phenotype and — more importantly — on how these genes work 
together, we developed a sensitive parallel technique for quanti- 
fying the CLS of yeast single- or double-gene deletions (Figure 1). 
We used fluorescent-protein labeling to track the apparent death 
rates of mutant and wild-type (WT) reference populations in 
stationary phase. Specifically, each mbced population suspended in 
stationary phase was automatically sampled at regular age 
intervals by inoculating to fresh medium, where competition 
assays were used to estimate the abundance of viable mutants 
relative to the WT at the time of sampling. An apparent relative 
lifespan, L, of each mutant strain was defined from the change in 
relative viability as a function of time in stationary phase. 
Importantly, we accounted for differences in growth rates during 
outgrowth (see model in Materials and Methods). By using fuUy- 
prototrophic strains in all screens, we avoided the strong effects of 
amino acid auxotrophy on stationary-phase survival [23,35]. In 
addition, we buffered the aging medium to identify the lifespan 
effects of genes and their interactions independentiy of acetic acid 
toxicity [36]. 

We evaluated the performance of our CLS method in several 
different ways. First, we estimated the technical and biological 
replicabHity of our L measurements. We observed that the results 
of two independent large-scale (720 deletion strains) experiments 



were highly correlated (Pearson's r=0.88; Figure SI). Dye-swap 
experiments in which 235 deletion strains were independendy 
tagged with both RFP and CFP fluorescent proteins also showed a 
good correlation (Pearson's r=0.79; Figure SI). In addition, we 
compared results of our competition-assay with an alternative 
high-throughput approach that is based on outgrowing individu- 
ally-cultured strains at different age intervals [20] (Figure S2). 
Competition-based monitoring of CLS resulted in considerably 
less variability in the measured phenotypes than experiments in 
which strains were individually aged and outgrown, underscoring 
the effectiveness of our experimental approach in terms of precision. 

Given that strain-strain interactions could obscure the results of 
a CLS screen based on competitive survival, we also evaluated the 
effects of using different reference strains. The relative lifespans of 
a set of 24 knockouts were mostly consistent across experiments in 
which different combinations of strains were used (/)<0.05, 
Spearman's correlation), indicating that the survival ranks were 
typically insensitive to the choice of reference strain (Figure S3). 

Taken together, these results indicate that our strategy possesses 
accuracy and robustness that enable the large-scale characteriza- 
tion of CLS at high resolution. 

A high-resolution genome-wide screen for lifespan 
factors 

To screen the genome for lifespan-determining genes and 
pathways, we used our CLS method to characterize 3,878 
nonessential-gene knockouts from the yeast deletion collection 
[37] (data set in Table SI). The distribution of lifespan effects 
peaked near neutrality {L = 1), with a long tail of short-lived 
knockouts and a smaller but substantial proportion of long-lived 
strains (Figure 2A). Given that our high-throughput screen was 
performed at low aeration and not under the conventional 
aeration conditions (see Materials and Methods), we characterized 
a subset of mutants at high aeration and confirmed that there was 
a general agreement in their long- and short-lived behaviors 
regardless of aeration (Figure S4). 

In genetic studies, mutations that increase a maximal lifespan 
are more likely to be informative about aging processes than 
mutations that reduce the life span [38], as such deleterious 
phenotypes can arise from a variety of reasons not necessarily 
related to aging. Indeed, we observed that short-lived knockouts 
are more likely to be affected in their rate of exponential growth 
than strains with wild-type CLS (Figure S5). Still, we note that at 
least 64% of the short-lived knockouts identified in our genome- 
wide screen were specific for CLS, that is, their deletion affected 
stationary-phase survival but not the rate of growth (Figure S5). 
Furthermore, the overlap of mutants affected in both stationary- 
phase survival and exponential-growth rate did not increase 
substantially when a subset of mutant strains were grown non- 
fermentatively under ethanol-glycerol medium (Figure S6), which 
indicates that respiratory integrity and its overlap with stationary- 
phase survival is largely taken into account by describing the 
mutant's fermento-respiratory growth capacity under standard 
glucose medium. 

We found that an important fraction of the non-essential yeast 
genome mediates CLS. Overall, 6.8% of single-gene knockouts 
had a long-lived phenotype, while 7.2% lived significantiy less than 
the wild type without showing major growth defects (5% FDR; see 
Figure 2a inset, genes in Table S2). To quantitatively describe the 
cellular processes and protein complexes that underlie CLS in 
yeast, we also identified functional groups of single deletions with a 
skewed distribution of lifespan effects (Figure 2B). Mutant strains 
with compromised mitochondrial function, de now purine biosyn- 
thesis, vacuole function, and autophagy typically resulted in 
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Cells in mixed cultures Populations die at different rates and A relative survival 

age in stationary phase the fraction of viable cells is estimated coefficient is obtained 

at different ages (slope) 



Figure 1 . A novel strategy enabling quantitative genetic analysis of yeast CLS. (A) Schematic showing an RFP-labeled mutant, xA, that lives 
shorter than the CFP-labeled WT reference. Strains are mixed in equal amounts and grown to saturation in 96-deep-well plates, where the mixed 
population ages. (B) At regular time intervals, outgrowth cultures are automatically inoculated from the same aging stationary-phase plate onto fresh 
medium and monitored for growth. The relative initial fraction of viable cells (magenta braces) is approximated from the change in relative 
fluorescence signal interpolated to a fixed time point (cyan vertical line) in exponential outgrowth. (C) The change in ratio of viable cells as a function 
of age is used to obtain an apparent relative survival coefficient, s, which is a direct estimate of the lifespan of xA relative to WT. 
doi:1 0.1 371/journal.pgen.1 0041 68.g001 



reduced lifespan phenotypes. This confirmed the main findings of 
recent genome-wide screens for yeast CLS factors [21-23]. We 
also note that genes involved in the regulation of replicative cell 
aging were enriched in CLS effects. 

An inspection of the list of short- and long-lived strains (Table 
S2) indicated that we were also able to confirm a number of 
previously-reported CLS genes. For instance, inactivation of the 
Leu3 transcription factor involved in amino acid biosynthesis 
increased CLS, as expected [39]. Lifespan extension was also 
observed for ald6A [40]. In another example, depletion of the 
apoptotic proteins Fisl and Aim 14 extended yeast CLS, in 
consonance with studies showing that chronological aging induces 
apoptosis in yeast [41,42]. This was also the case of the 
mitochondrial ubiquinone oxidoreductase: disruption of Ndil 
decreases ROS production and elongates the CLS of yeast [43]. 
Moreover, strains impaired in Msn2, Glc8, and Glc3 had reduced 
lifespan, which was in agreement to studies showing that decreased 
Tor activity activates stress-response transcription factors and 
genes for glycogen accumulation [19]. 

Importantly, our screen allowed us to uncover biological 
processes and protein complexes that had not been previously 
associated to the CLS phenotype. Noteworthy, depletion of a 



number of chromatin-modification proteins resulted in extended 
lifespan (Figure 2B). Such novel chronological aging factors 
included the Swrl -complex proteins Arp6, Swc3, and Swrl. 
The SAGA-complex subunits Chdl and Sgfl 1 were also, to our 
knowledge, previously unrecognized pro-aging proteins, which 
add to proteins of the SAGA-like complex as a histone- 
modifying factors underlying yeast lifespan [44]. On the other 
hand, disruption of the Rpd3L histone-deacetylation complex 
subunits Cti6, Depl, Rxt2, Rxt3, and Umel resulted in short 
lifespan, without altering the cell's growth capacity. Likewise, 
the Sdcl, Sppl, and Swdl subunits of the Setl/COMPASS 
histone-methylation complex were novel yeast CLS factors (see 
Table S2). We also identified a number of novel lifespan 
phenotypes associated to genes with roles on sterol and 
sphingolipid metabolism and homeostasis (Figure 2B; Table 
S2). Depletion of the ergosterol-biosynthetic genes Erg5, Erg6, 
and Hmgl extended CLS, while disruption of the fatty-acid 
elongase Surl — involved in sphingolipid biosynthesis — reduced 
CLS without affecting growth. The ARVl -encoded protein, 
which mediates sterol and shingolipid homeostasis [45-48], was 
also a novel lifespan factor identified in our genome-wide CLS 
screen. 
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Figure 2. High-resolution genome-wide screening of yeast CLS. (A) Cumulative frequencies of the relative lifespan of most viable knockout 
strains (black circles, /i = 3,878) and replicates of the WT (gray crosses, n = 414). Colored symbols correspond to selected knockouts (right panel). Inset 
shows the fraction of significant short- (magenta) and long-lived phenotypes (cyan) (5% FDR). Short-CLS phenotypes with reduced growth capacity 
(G<0.95) were filtered-out from subsequent analyses (see Figure S5). (B) CLS-specific gene groups and pathways; plot shows the median relative 
lifespan (L) of genes annotated in the same GO term, which was compared in each case to the median lifespan of the entire data set. Labels at the 
vertical axis show representative enriched GO categories (p<0.05, Wilcoxon rank sum test). The radius of each circle indicates the number of 
measured gene knockouts within the GO group, the minimum being five (de novo IMP biosynthesis) and the maximum 39 (nucleolus); colors were 
arbitrarily assigned to similar functional categories. Data set includes 3,286 strains that have no marked defect in growth (G>0.95). 
doi:1 0.1 371 /journal.pgen.1 0041 68.g002 
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Confirming chronological lifespan phenotypes and 

estimating extents of adaptive regrowth 

To further validate the resuhs of our competition-based 
genome-wide screen, we assayed 14 of the longest-lived and seven 
of the shortest-Uved knockout strains using the standard CLS 
method with non-buffered medium [49] (see Materials and 
Methods). Knockout strains Arps4b, ATps24a, ArpslOb, /S.rpslSh, 
Aku3, Arcyl, and Afl/(/6'were confirmed as long-lived (Figure 3A- 
G), while Alscl, Aglkl, Alsc2, Astfl, and Auhpl 4 were confirmed as 
short-Kved (Figure 3H-L). The CLS of strains deleted for ADE5,7, 
AQRl, CrB5, ERG6, JVEWl, RRP8, SPPl, TAKl, and rPL062W 
was not distinguishable from the WT when retested (data not 
shown). Overall, seven out of 14 long-lived strains and five out of 
seven short-lived knockouts were confirmed using the standard 
CLS method. 

Subpopulations in aging stationary phase cultures may undergo 
adaptive regrowth, a phenomenon whereby non-dividing cells 
acquire mutations that allow them to re-enter the cell cycle and 

take advantage of the nutrients released by dead cells [.50]. If the 
rate of adaptive regrowth of a given mutant strain was different 
than the wild-type, such mutant would register an artificially high 
or low abundance ratio as a function of age. We therefore 
estimated the extent at which adaptive regrowth gave rise to false- 
positive lifespan factors in our screen (Figure S7). To do so, we 
focused on nearly all mutant strains scored as short- or long-lived 
(505 out of 511 knockouts; see Table S2) and monitored their 
chronological aging in water to prevent the occurrence of adaptive 
regrowth [49]. We observed that the phenotypes from these two 
aging conditions were consistent with each other; the distributions 
of lifespan effects in water were markedly different between 
mutants scored as short-lived or king-lived in buffered medium 
{p<lO WUcoxon rank sum test). In other words, strains that 
aged fast or slow in buffered medium typically did so in water as 
well. For instance, the CLS of only one of the longest-lived and 
one of the shortest-lived strains in buffered medium were not 
significantly diflFc-rcnt from the WT CLS in water (see Table 1). 
Overall, we were abli; to confirm 61% short-lived and 48% long- 
lived phenotypes, which sets such minimal fractions of true- 
positive lifespan factors in our original genome-wide screen in 
terms of accounting for adaptive regrowth. 

The Swrl complex is a novel aging factor that mediates 
lifespan extension by dietary restriction 

Swrl is a multisubunit chromatin-remodeling complex that 
contains the ^'IfRZ-encoded Swi2/Snf2-related ATPase and is 
required for the incorporation of the histone variant H2A.Z into 
chromatin [51]. To further investigate the mechanisms by which 
activity of the Swrl complex promotes chronological aging, we 
quantified the lifespan effects of Swrl mutants under glutamine 
(nitrogen-rich) or GABA (nitrogen-poor) conditions (Figure 4). By 
doing so, we asked whether lifespan extension by dietary 
restriction depends on the reduced activity of this novel aging 
factor. Interestingly, the Swrl -complex gene knockouts hved 
relati\ely longer than the WT in glutamine, but not under 
nitrogen-limited condition, with the exception ot arplA which had 
a significant effect in both GABA and glutamine (Figure 4). As a 
control, we analyzed the effects Tori and Rim 15, which are 
known to be key players of the pathway that extends longevity by 
nutrient limitation [13]. As expected, the tori A and rimlSA strains 
showed, respectively, increased and decreased relative longevity 
preferentially under nitrogen-rich conditions. These results suggest 
that lifespan extension by dietary restriction depends, at least in 
part, on the action of the Swrl histone-exchange complex. 



A functional link between Arvl and autophagy underlies 

yeast lifespan 

To gain deeper knowledge on the lifespan role of Arvl, a 
protein involved in lipid homeostasis, we asked whether this novel 
CLS factor modulates yeast lifespan together with or indepen- 
dentiy from autophagy, which is a conserved modulator of lifespan 
in yeast and multicellular eukarv'otes [14,19,52] and is intimately 
linked to lipid metabolism [53]. We generated a collection of 
double knockouts combining arvlA to 26 deletions of core- 
autophagy or CVT-pathway genes. Strikingly, the relative lifespan 
of many double knockouts was significantiy increased when 
compared to the expected from the independent combination of 
the deleterious lifespan phenotypes of each of the corresponding 
single knockouts (Figure 5A). This result suggested that the anti- 
aging effect of autophagy depends to some extent on the action of 
Arvl. 

To confirm the observed genetic association of Arvl with the 
autophagy machinery, we monitored different reporters of 
autophagic activity in wild-type and Arvl -disrupted cells. We 
found that awlA cells accumulate less Atg8 — a marker of 
autophagic bodies — in the vacuole, partially resembhng the 
molecular phenotype of atglA cells in which autophagy activity 
is completely blocked (Figure 5B). In addition, we examined the 
degradation of Pgkl-GFP, a low-turnover rate cytosolic fusion 
protein that is rapidly degraded by non-selective autophagy upon 
nitrogen starvation [.54]. Also in agreement with the genetic 
interactions scored between autophagy genes zndARVl, we found 
that the autophagy-directed appearance of processed OFF protein 
was clearly diminished in arvl A cells (Figure 5C). Taken together, 
these results showed that autophagic activity was compromised in 
arvl A deletion strains and confirmed a novel connertion of Arvl 
with the autophagy pathway, as anticipated by means of genetic- 
interaction analysis. 

Deriving genetic interactions systematically from 
quantitative lifespan phenotypes 

The combination of high throughput and sensitivity of our 
phenot^ping approach not only provides an accurate catalogue of 
genes and pathways underlying CLS, but also lends itself for the 
systematic construction of lifespan-based quantitative genetic 
interaction (GI) maps. By describing the ways in which mutations 
affect each other's phenotypes, GI analyses in yeast have been very 
successful in elucidating genetic pathways [27-31,55]. Since cell- 
proliferation phenotypes are typically used to generate such maps, 
and given that the genes that promote growth are not neci;ssarily 
important for survival, it was expected that CLS-based GI data 
would uncover novel pathways, with particular relevance to aging. 
Indeed, we noted that 64% of the genetic factors identified in our 
genome-wide screen were specific for stationary-phase survival, 
that is, their deletion affected CLS but not the rate of exponential 
colony growth (Figure S5). 

For the systematic study of lifespan-epistatic interactions we 
focused on the autophagy pathway, which drives the degradation 
and recycling of unnecessary or dysfunctional cellular components 
[,56-58]. While autophagy has consistentiy been shown to 
modulate CLS, the way in which autophagy genes interact 
phenotypicaUy with one another and with genes of related 
processes has not been systematically analyzed. We surveyed 
interactions for almost all core-autophagy genes, along with 
several signaling and accessory factors (39 genes and all possible 
combinations). The single knockout of most genes examined had 
at least a small but significant deleterious effect on CLS, and the 
measured CLS phenotypes of two marker-swap experiments 
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Figure 3. Longevity plienotypes were confirmed using tKie standard CLS method. Survival curves under non-buffered SC medium for WT 
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doi:1 0.1 371 /journal.pgen.1 0041 68.g003 



showed a good correlation (Figure S8). Importantly, to verify the 
CLS specificity of the genes tested, we confirmed that exponential 
growth capacity of their single knockouts remained largely 
unafiected (Figure S9). 

For each gene pair, lifespan epistatic interactions were defined 
as significant deviations from a null multiplicative expectation 
based on the CLS phenotype of the two single knockouts. 
Neutrality was defined when the lifespan plienotypes of the single 



knockouts combined multiplicatively, suggesting that the two genes 
act independently from one another (Figure 6A, left panel). 
Departing from this neutral expectation, GIs were classified as 
either negative (aggravating) or positive (alleviating) [59]: Negative 
epistasis was scored when the double mutant exhibited a 
significantly shorter CLS phenotype than one would expect from 
the CLS of the two single knockouts (Figure 6A, center panel). 
Such aggravating phenotype is expected from the deletion of pairs 
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Table 1. Chronological-lifespan effects of strains aged in SDC medium and In water. 



SDC buffered medium 



Water 



Strain 


Relative lifespan, L' 


p-value"" 


Relative lifespan, L' 


p-value"" 


WT 1 1 


Long-lived strains 


rpsAbA 


1.033±0.003 


<1C"^ 


1.009 ±0.0 10 


0.105 


rps24aA 


1.035±0.007 


<1C"^ 


1.014±0.0C9 


0.045 


rpsWbA 


1.042 ±0.004 


<10"^ 


1.020±0009 


0.006 


rpslSbA 


1.035±0.003 


<10"^ 


1.014±0.014 


0.043 


leuSA 


1.01 2 ±0.004 


0.003 


1.01 5 ±0.006 


0.032 


rcylA 


1.062 ±0.006 


<1C"^ 


1.051 ±0.005 


<10"^ 


ald6A 


1.1 09 ±0.005 


<1C"^ 


1.063 ±0.005 


<10"^ 


Short-lived strains 


IsclA 


0.964±0.003 


<10"^ 


0.975 ±0.009 


<10"^ 


glklA 


0.957±0.003 


<10"^ 


0.999 ±0.0 15 


0.444 


bc2A 


0.964±0.002 


<1C"^ 


0.984±0.006 


0.022 


stflA 


0.959±0.003 


<10"^ 


0.973±0.007 


<10"^ 


ubp 14A 


0.963 ±0.004 


<10"^ 


0.971 ±0.010 


<10"^ 



'^Relative lifespan is obtained from the survival coefficient (l+s), expressed as the change in viability per day relative to the WT reference, ± error in the linear fit. 
''p-value is from the Z-score of 414 or 65 WT replicates in SDC medium or in water, respectively. 
doi:1 0.1 371 /journal.pgen.1 0041 68.t001 



of genes that act in lifespan pathways that partially compensate for 
each other's loss. Conversely, in a positive epistatic interaction, the 
double mutant exhibited an extended CLS relative to the neutral 
expectation (Figure 6A, right panel). In general, this phenotypic 
response takes place when both deleted genes are part of the same 
linear lifespan-determrning pathway or protein complex. 

Overall, we scored a significant GI for 20 1 out of 72 1 gene pairs 
tested (3.2% negative and 24.7% positive, 95% CI, see data in 
Table S3). In agreement with the common assumption that most 
gene pairs do not interact, we found that the lifespan-epistasis 
spectrum peaked close to neutrality, though skewed to positive 
epistasis (gray histogram in Figure 6B). Many of the positive 
interactions in our data set corresponded to gene parrs within the 
core-autophagy machinery (purple histogram outline in Figure 6B). 
For this set, we scored 55.8% significant positive interactions and 
no negatively interacting gene pair. The distribution of epistatic 
effects remained largely unaffected when interactions were defined 
from an additive rather than multiplicative neutral expectation 
model (Figure SIO). 

We also observed that the distribution of epistasis among gene 
pairs that code for proteins that interact physically had a higher 
median s value than that of gene pairs whose protein products do 
not interact (Figure 6C). In summary, the observed prevalence of 
positive epistasis within our data set was consistent with the fact 
that many autophagy proteins operate in concert as a non-essential 
functional unit that promotes CLS. 

A lifespan-epistasis map for autopliagy 

Based on the correlation of the quantitative epistatic-interaction 
pattern derived from CLS phenotypes of each gene, we 
constructed a high-density GI map for autophagy (Figure 7A). 
Using purely phenotypic information, this unsupervised map 
recapitulated known pathways and complexes and shed hght on 
novel gene associations. With few exceptions (Atg5, Atgl2, and 
Vps30), all core autophagy proteins were grouped together (core- 
ATG module. Figure 7A). Proteins with a major role in the CVT 



pathway (Atgl9, Atg23, and Atg27) were clustered together with 
other autophagy-accessory proteins and apart from the core-ATG 
module (CVT module. Figure 7A). 

We also used gene-clustering information along with gene- 
interaction significance to generate a lifespan-epistasis network 
(Figure 7B). We observed a high enrichment of positive epistatic 
interactions among genes within the core-ATG module, again 
consistent with the fact that core-autophagy proteins act as a single 
functional unit (p<lO~*, Fisher's exact test). Furthermore, 
directionality was assigned to positive epistatic interactions when 
one of the genes had a larger single-knockout effect than its 
interaction partner [30]. For instance, we scored positive 
interactions targeted to ATG18; single-deletion of this gene had 
a smaller quantitative effect on CLS when compared to the typical 
core-autophagy gene. This observation was consistent with the fact 
that Atgl8 is only required for reaching full autophagic activity 
and most autophagy proteins localize to the pre-autophagosome in 
the absence of Atgl8, but not in the absence of other core- 
autophagy proteins [60] . 

Our lifespan-epistasis map also recapitulated a group of proteins 
related to phosphatidylinositol phosphate metabohsm and regula- 
tion (Figure 7B), including Atg21, Tax4, Snx4, Hsv2, and Tepl. 
The association of Tep 1 within the PIP module provided genetic 
evidence for its putative role as a phosphornositide phosphatase. In 
cancer cells, the TEPl homologue PTEM controls autophagy by 
downregulating the PI3KclassI/PKB pathway [61]. The syntaxin- 
like t-SNARE Tlg2 protein was also part of this PIP module. 
Recent evidence suggests that vesicle-tethering proteins and some 
SNARE proteins are required for autophagosome biogenesis and 
normal anterograde transport of Atg9 [62]. This result provided 
further genetic evidence for the role of Tlg2 as a SNARE protein 
which is essential for autophagy and suggested that this functional 
association is needed for normal lifespan. 

Three signaling factors (Pho80, Snfl, and Regl) were grouped 
together in our lifespan-epistasis map (Figure 7A,B). For PHO80, 
positive epistatic interactions were all directed from this signaling 
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Figure 4. Relative lifespan extension by inactivation of Swr1 is partially abolished under dietary restriction. Bars sl-iow tl-ie relative 
lifespan effects {L—^) for mutant strains aged in competition against the WT in glutamine (rich nitrogen, blue bars) or GABA (poor nitrogen, orange 
bars) aging medium. Error bars indicate the standard error from at least three independent experiments. Box plot (inset) shows the absolute half life 
measurements of individual WT cultures aged under glutamine or GABA (n = 28). 
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gene to autophagy genes, in agreement with the regulatory role of 
Pho80-Pho85 upstream of autophagy [63] and suggesting that this 
cyclin regulates other downstream pathways that are important for 
yeast CLS. As opposed to PHO80, SNFl interacted negatively with 
the core-ATG and PIP modules, which may reflect that Snfl 
independently promotes autophagy and inositol biosynthesis 
[64,65]. 

Taken together, these results show that quantitative lifespan- 
epistasis maps represent a powerful systematic strategy to estabhsh 
functional associations within and between aging pathways. 

Discussion 

We have used a novel functional genomics platform to generate 
a high-resolution compendium of genes that affect yeast station- 
ary-phase survival and to construct the first systematic GI map 
based on lifespan phenotypes. 

Our genome-wide screen indicated that a substantial fraction 
of the non-essential yeast genome (14%) regulates stationary- 
phase survival independently of general growth defects. Among 
novel CLS factors detected in our screen were different 
proteins of the Swrl complex, which is required for the ATP- 
dependent recruitment of histone variant Htzl — the yeast 
homologue of mammalian H2A.Z — and promotes expression 
near silent heterochromatin [51]. It is known that chromatin 
structure and dynamics change during the aging process [66]. 
To further explore the specific role of the Swrl complex as 
a regulator of CLS in yeast, we characterized a number of 



Swrl -complex knockouts under different nitrogen conditions. 
In doing so, we provided evidence suggesting that lifespan 
extension by dietary restriction depends on the activity of the 
Swrl complex. 

We also showed that Arvl — an endoplasmic reticulum protein 
conserved from yeast to mammals — and autophagy promote 
chronological longevity in concert. Such genetic association of 
ARVl was confirmed by monitoring different autophagy reporters 
in Arvl -depleted cells. This result adds to a growing body of 
evidence suggesting that autophagy is intimately connected to lipid 
metabolism: While autophagy regulates lipid metabolism [53], 
sphingolipids and ceramides play important roles in autophagy 
itself, both at the level of regulation and as membrane constituents 
of autophagic vesicles [67,68]. In the absence of Arvl, yeast cells 
up-regulate the unfolded protein response [69] and display a 
variety of lipid-related phenotypes that include aberrant sterol and 
sphingolipid metabolism and loss of plasma-membrane asymmetry 
[45-48]. In humans, Arvl is required for normal cholesterol and 
bile acid homeostasis [70] . In regard to our finding that autophagy 
activity depends on Arvl, we speculate that this protein could be 
specifically required for proper targeting of Trs85 to the pre- 
autophagosomal structure, as it has been recently shown that Arvl 
is involved in the targeting of tail-anchored proteins to the 
endoplasmic reticulum [48]. Alternatively, the general lipid- 
homeostasis role of Arv 1 could be essential for the correct activity 
of the autophagy machinery. In any case, our results indicate that 
the functional link between Arvl and autophagy has an impact on 
cellular lifespan. 
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Figure 5. Functional link between autophagy and Arv1 underlies CLS. (A) Genetic interaction analysis of ARV1. The relative lifespan of 26 
arvlA, xA double knockouts is plotted against the expected multiplicative combination of the corresponding pairs of single knockouts. Dotted lines 
at the horizontal and vertical axes show the effect of the arvlA single knockout. Green circles indicate significant deviations from the neutral 
expectation (>95% CI), which are shown as green edges connecting ARVl to the interacting gene (right). (B) Micrographs of WT and gene-knockout 
strains expressing the GFP-Atg8 fusion protein under basal conditions (SD, left) and after 7 h of nitrogen starvation (SD-N, right). Vacuoles were 
stained with FIVI4-64; DIC images show the position of yeast cells. (C) Western blot analysis of nitrogen-starved strains expressing the fusion protein 
Pgkl-GFP: WT and cells depleted for Arvl, the core-autophagy Atgl protein, or the accessory Atg21 protein. Time zero corresponds to mid-log 
cultures before transfer to SD-N medium where cells grew for 2, 5, or 7 h before lysis. 
doi:1 0.1 371/journal.pgen.1 0041 68.g005 



Our systematic lifespan-epistasis network successfully recapitu- 
lated functional associations of the autophagic machinery and 
related biological processes. While most core-autophagy genes 
clustered together as expected, few exceptions included ATG5 and 
ATG12, which were grouped with genes of selective autophagy 



and of the CVT pathway. A specific role of ^TG5 and ATG12 has 
been reported for mitophagy and lifespan regulation in human 
umbilical vein endothelial cells: Mitochondrial damage induces 
upregulation of these genes together with LC3B/ATG8, which in 
turn increases the replicative cellular lifespan [71]. Furthermore, 
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Figure 6. Lifespan-epistatic interactions recapitulate functional and physical associations. (A) Plots show the measured relative lifespan 
of representative gene pairs (single and double knockouts) with neutral (left panel), negative (center panel), and positive epistasis (right panel). 
Dashed lines indicate the multiplicative expectation in each case, that is, L^a • i-yA- Error bars indicate the 95% CI. (B) Spectrum of epistasis measured 
for all gene pairs (gray histograms, n = 721) and for core-autophagy gene pairs (purple histogram outline, n = 75); (C) Comparative cumulative 
distributions of 8 for gene pairs coding for proteins that interact (orange line, n = 31) or do not interact physically (black line, n = 690); the p-value is 
from the Wilcoxon rank-sum test. 
doi:1 0.1 371/journal.pgen.1 0041 68.g006 



our epistasis network established that Tepl and Tlg2 work in 
concert with the pliosphatidylinositol-phosphate metabohsm 
pathway to modulate lifespan. While it has been shown tliat 
Tepl is not a replicative lifespan factor in yeast [72], our data 
indicated that Tepl activity promotes CLS, in line with PTEN- 
overexpressing mice witli increased longevity [73]. The CLS 
phenotype of yeast Tep 1 and its association to other phospliati- 
dylinositol proteins provide a simple model system to elucidate the 
molecular mechanisms underlying the role of its mammalian 
homologue PTEN in embryonic survival and tumor suppression 
[74]. 

Compared to previous genome-wide screens in which 6% [19] 
or 31% [21] of longevity phenotypes were confirmed when 
retested, our approach performed modestly better by correctly 
predicting 50% (seven out of 14) knockout strains with extended 
lifespan. There are a number of possible reasons why not all of the 
mutants retested with the standard CLS assay were successfully 
confirmed. Our relative lifespan measurements are based on a 
simple deterministic model that assumes constant exponential 
death and outgrowth rates as the population ages. For instance, 
cell heterogeneity of stationary-phase cultures [75] was not 
possibly taken into account. Moreover, results of our control 
screens suggest that difiFerences in aeration or the effects of 
adaptive regrowth explain part of the false positives that were 
obtained. Importantly, once a set of mutants have been confirmed, 
our approach can be readily used to leverage the chronological 



aging paradigm to generate aging networks, as we have shown for 
autophagy. Our pure phenotypic characterizations proved pow- 
erful for establishing functional relationships among such genes 
that had remained hidden to large-scale genetic analysis. Despite 
the fact that our epistasis data set is limited in its scale and is a priori 
enriched for genes that share functional relationships, we 
anticipate that the strategy herein described will open the door 
to describe the spectrum of lifespan-epistatic interactions in a 
comprehensive and unbiased manner. 

With a growing number of lifespan-related genes, a systems view 
of aging is much needed to grant a deeper mechanistic 
understanding of cellular longevity [7,76]. Our quantitative anal- 
ysis of longevity factors and their genetic interactions in yeast is an 
important step to this end. 

Materials and Methods 

Strains and media 

Fluorescent SGA starter strains YEGOl-CFP and YEGOl-RFP 
[AUTa. PDCl-XFP-CaURA3MX4 canlA::STE2pr-SpHIS5 lyplA 
hi.s3Al uraSAO LEU2) were generated by direct PCR-based 
tagging of the S288C-derivative SGA starter strain Y8205 [77] 
with the pBSlO Gem\ea-n-hphMX4 (CFP) or pBS35 mCherry- 
hphMX4 (RFP) cassettes (The Yeast Resource Center) at the PDCl 
locus and subsequent direct replacement of the hphMX4 cassette 
with CaURA3MX4; strong constitutive expression of fluorescent 



PLOS Genetics | www.plosgenetics.org 



10 



February 2014 | Volume 10 | Issue 2 | e1004168 



Systematic Genetic Analysis of Yeast Longevity 



I 



0.1 



-0.1 




core ATG 



CVT 



signaling 



PIP 



B 



9 Core ATG 

CVT pathway / accessory ATG 
PIP metabolism / regulation 
ATG signaling 



s<0 

S>0, Lx = Ly 
£>0, Lx > Ly 




Figure 7. Quantitative lifespan-epistasis analysis reveals functional associations within and between pathways underlying CLS. (A) 

Hierarchical-cluster analysis of lifespan-epistasis data. Each column/row represents the interaction pattern for a specific gene; the color scale indicates 
the epistasis interaction strength, e. Dendrogram clusters are highlighted in colored bars (right). (B) Network representation of gene clusters and their 
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Arrows indicate cases of significant CLS-differences between the two single knockouts, with directionality from the largest to the smallest single 
knockout effect. 
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proteins during exponential growtli was acliieved by carboxyl- 
terminal fusion to the Pdc 1 protein. Deletion strains were from the 
yeast deletion collection (Open Biosystems) in the S288C- 
derivative BY4741 background [A'lATa xxxA::kanMX4 hisSAl 



uraSAO leu2A0 met 15 AO). The GFP-ATG8 fusion was generated 
by PCR-based tagging with the natMX4-CUPlpr-yeGFP cassette 
from pYM-N4; the PGKl -GFPimion was obtained from the GFP- 
tagged yeast strain collection (Invitrogen). Genes of interest were 
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knocked-out from GFP fusion strains by direct gene replacement 
with the kanMX4 cassette. 

The following growth media were used: (1) Aging medium: 
synthetic complete (SC) medium with 2% glucose, 0.2% amino 
acid supplement mbc as defined (Cold Spring Harbor Laboratory 
Manual 2005), and buffered to pH 6.0 with 17.9 mM citric acid 
and 64.2 mM dibasic sodium phosphate. When noted, 25 mM 
glutamine or 25 mM gamma-aminobutyric acid (GABA) instead 
of ammonium sulfate and 0.07% of amino acid supplemented mix 
were used. (2) Low-fluorescence medium (YNB-lf) [78]. (3) 
Nitrogen-star\'ation medium (SD-N): 0.17% yeast nitrogen base 
without amino acids and ammonium sulfate and with 2% glucose. 
(4) Yeast-extract peptone with 2% dextrose (YPD) or 2% ethanol, 
2% glycerol (YPEG) medium. 

Generation of yeast libraries 

Collections of prototrophic yeast strains were generated by 
synthetic genetic array (SG.4) methodology [77]. Colony arrays 
were transferred manually with a 384-head pin tool (V&P 
Scientific, VP384F); antibiotic concentrations used for selection 
were 200 |tg/ml G418 (Invitrogen) and 100 ^lg/ml clonNAT 
(Werner Bio Agents). 

To generate the collection of RFP-tagged strains used for the 
genome-wide CLS screen, the YEGOl-RFP starter strain was 
mated to an array of 4,844 viable single-deletion strains, followed 
by diploid selection, sporulation, and three rounds of haploid 
selection [HI^ for MATa mating type, URA^ for fluorescence 
marker, and G418^ for knockout selection). The resulting 
prototrophic strains were MAT^ PDCl-RFP-CaURA3MX4 
canlK::STE2pr~SpHIS5 lyplA uraSAO his3M LEU2 METIS 
xxxA::kanMX4. 

For the double-knockout collection of autophagy-related genes, 
41 query strains were generated by PCR-based gene replacement 
of the YEGOl-RFP background using the natMX4 module. This 
query collection was mated to an array of 44 deletion strains. 
Double-knockout haploid strains were selected as described above, 
plus clonNAT*" for the second knockout. The query collection 
included two neutral-insertion AoA controls, while four hisSA 
controls were present in the array collection. All combinations of a 
knockout strain with a neutral insertion (up to six replicates) define 
the "single knockouts", while the combination of two neutral 
markers define the "WT strain". Likewise, the CFP-labeled 
reference strain was obtained from the YEGOl-CFP starter strain 
by the same procedure of combining two neutral markers. All 
query and array strains were PCR-verified for correct marker 
insertion and disruption of the WT allele (data not shown). 

To obtain ARVl double knockouts, awlA::natA4X4 and hoA 
neutral-insertion query strains (YEGOl-RFP background) were 
mated in triplicate to an array of 29 autophagy-gene deletion 
strains, including the arvlA::kanMX4 and the his3A::kanMX4 
neutral-insertion. Double knockouts were selected as described 
above. 

Automated CLS assay: Experimental setup, robotic 
integration, and data acquisition 

Saturated cultures of RFP-labeled mutant (single or double 
knockout) or CFP-labeled reference strains grown in SC medium 
were mixed in a 1:1 ratio to a final volume of 150 |J,1 in 96-well 
plates (Corning 3585) and pinned with a manual replicator (V&P 
Scientific, VP 407) onto 700 |il of fresh aging medium in deep-well 
plates (Nunc 26025 1) with disposable plastic covers. All subsequent 
steps were carried out in an automated robotic station (Tecan 
Freedom EVO200) that integrates a plate carrousel (Liconic 
STXllO), a high-performance multilabel plate reader (Tecan 



Infinite Ml 000), a 96-channel pipetting head, an orbital plate 
shaker, and a robotic manipulator arm, all contained in an isolated 
environmental room. Mixed stationary-phase cultures were 
maintained at constant temperature (30°C) and relative humidity 
(70%) without shaking; cultures were fully resuspended every 24 h 
by automated pipetting. 

Five days after inoculation (age zero, as previously defined [50]) 
and every 2-3 days for up to 23 days from this point, stationary- 
phase cultures were transferred to the hquid-handling deck for 
outgrowth: 5 )J.l aliquots were inoculated into 150 |J,1 of fresh low- 
fluoresccnc(; medium in 96-well plates. Such outgrowth cultures 
were kept at 30°C without shaking and monitored by the following 
procedure: Every 3 h (15 h, five data points), outgrowth cultures 
were vigorously shaken and transferred to the plate reader to 
collect data for both fluorescence channels (CFPraw: Ex 433 nm/ 
5 nm and Em 475 nm/5 nm; RFPraw: Ex 587 nm/5 nm and Em 
610 nm/5 nm) and optical density at 600 m {ODggg). Background 
fluorescence signal as a function of ODggg was routinely collected 
from RFP-only (CFPbg) and CFP-only {RFPbg) control cultures. 
Fluorescent signal in the outgrowth culture was defined as: 
RFP=RFPraw-RFPbg and CFP= CFPram-CFPbg. 

Data analysis 

The log of the ratio of AFP to GFP signal {InR/Q was obtained 
and a single value for each outgrowth culture at age T was 
calculated by interpolating InR/ C to a fixed point in time, t* (see 
Model below). For all experiments, t* was fixed at 10 h. 
Importandy, results were insensitive to the choice of t* for 
outgrowth data collected during late exponential growth (8-12 h 
post inoculation, data not shown). An apparent survival coefficient, 
s, and its standard error, CT^, were obtained from the slope of the 
linear fit of InR/Cf to T (robust regression, Madab). From an 
empirically determined dynamic range (Figure Sll), we defined 
conservative low and up InR/Ct* thresholds of —3.5 and 2.5, 
respectively. Only fits with at least three data points were 
considered to obtain a value of s. Hence, mutant strains with 
low fluorescence signal or strains that grew too slow in co-culture 
(less than 5% of the tagged deletion collection) were excluded from 
our analyses. 

All screens included a large number of WT samples scattered 
through the strain arrays (typically >300 wells); these were used to 
account both for survival differences between WT and CFP- 
labeled reference strain and potential batch effects. By definition, 
the survival coefficient of the RFP-labeled WT strain, Sv,t, equals 
zero. Hence, data was normalized by subtracting the median value 
of all WT control samples, Swt, to all raw s values in each batch (see 
example in Figure SI 2). The mutant's apparent relative lifespan 
was expressed as the decrease in mutant population relative to 
WT, diat is, L^ = \+s^. 

Model for outgrowth and stationary phase populations 
in co-culture 

In our experimental setup, outgrowth from stationary phase is 
used to estimate the relative numbers of viable cells across an 
experiment. In practice, this quantity is obtained after several 
hours of inoculation into fresh medium. In what follows, we show 
that measuring the ratios of two exponentialh'-gTo\ving popula- 
tions at a fixed time-point in hours, f , across different ages of the 
culture in stationary phase, T, decouples potential differences in 
growth rates from differences in the quantity of interest, that is, the 
rates of death. In other words, the decrease in relative abundance 
of a slow-growing strain will remain constant and will thus not 
affect the estimation of a change in relative viability (slope), 
provided that such relative abundance is measured in each 
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outgrowth iteration at the same time after inoculation onto fresh 
medium. 

Let Nx and JV„, be the sizes of two populations, a mutant (x) and 
a wild-type reference [wt), growing exponentially in co-culture at 
rates & and g„{. 



-gxt 



Equation 2 with the relationship between the death rate and the 
half life of a population, A = — and obtaining an estimate of 
the mutant's half life: 



ln(2) 



ln(2) 



Therefore, the relative number of viabk; cclh from an outgrowth 
curve inoculated at age T (measured in days) in stationary phase 
and measured at a fixed time point after inoculation (f, measured 
in hours), is given by: 



in 



(r,f)=ln 



X,tQ 



(1) 



What we get is our measure of interest — the relative number of 
viable cells at age T — plus a constant term that depends on a 
potential difference in growth rates and which can be ignored. 
This convenient decoupling would not be possible if the 
measurement was taken, for instance, at a fixed optical density 
of the culture. We note that, experimentally, population sizes were 
large (>10^) and all data was collected during exponential growth 
(see above). 

To model aging populations, we consider the simplest model in 
which mutant (x) and reference [wt) populations in co-culture die at 
constant exponential rates [79], and r^, resulting in population 
sizes jV^ and that decay exponentially as a function of age, T: 



-<-xT 



To measure the relative number of viable cells of each population, 
we take the logarithm of the population ratios as a function of age 
and obtain the equation for a line, with slope equal to the 
difference in the death rates of the two populations: 



In 



(r)=ln 



+ {rw 



(2) 



where .t is the apparent survival coefficient and In [l^J (-^) 

quantity that is obtained experimentally by estimating the number 
of viable cells at different ages of the culture (Equation 1). Since 
the death rate of the mutant, r„ is positive, the death rate of the 
reference strain, r^,, sets an upper bound for the value of ,s at s^r„t. 
We note that, experimentally, only a small number (16 out of 
3,878) of CLS measurements exceeded this theoretical upper 
bound, presumably due to the higher error in the estimate of 
extreme positive values (>t = NaN; Table SI). 

Where needed, we estimated an absolute apparent lifespan 
value for the mutant from its relative s value by combining 



where X^i is the measured half life of the RFP-labeled WT strain. 

Genome-wide screen and determination of lifespan 
factors 

A total of 4,340 mutant strains from the yeast deletion 
collection were successfully tagged; 4,050 reached ODggo and 
RFP signal above an arbitrary threshold. This collection was 
rearranged in 47 96-well plates which included 414 WT-RFP 
strains as controls and 48 additional wells with fresh medium or 
individual cultures of each WT-RFP or CFP-reference to account 
for background signal. The RFP-labeled deletion collection was 
screened against a WT CFP-reference strain with a half-life of 
21.7 days. After data acquisition and analysis, a relative survival 
coefficient (corrected s) was successfully calculated for 3,878 
strains (data in Table SI). 

Significandy short- and long-lived strains were scored by 
assigning a Z-score (^ to each knockout obtained from its L 
value and the standard deviation of 414 L measurements for the 
WT (normally distributed, p = 0.48, x^ test). Two-tailed /(-values 
were obtained from J^to estimate the false discovery rate (FDR) for 
multiple hypothesis testing; significant phenotypes were assigned 
using a ^<0.05 cutoff. 

Functional enrichment analysis 

The median relative lifespan (L) of genes from c-ach Gene 
Ontology (GO) term was compared to the median relative lifespan 
of the entire data set (i,= 1). GO terms with an L value which was 
significandy different from that of the genome are shown in 
Figure 2B (/)<0.05, WUcoxon rank sum test). GO annotations 
were downloaded from the Yeast Genome Database (April 2013); 
only those GO terms with at least five and no more than 60 genes 
and genes with a successful L measurement and relative growth at 
O0.9 were considered for analysis. 

Double knockout screen, definition of epistasis, and 
hierarchical clustering 

The final collection of 1,845 double-gene knockouts (41 query 

strains X 45 array strains) included up to six single-knockout 
replicates per gene of interest and two marker-swap replicates per 
double knockout (see correlation in Figure S8). Such indepen- 
dently-measured replicate values of L±a were combined to obtain 

an error-weighted average relative lifespan, L = 2__/ -i ^2^' 

with standard error ctl = < L^:^ 

V E,=iff, 

For every double deletion of genes X and T, epistasis was 
defined from a multiplicative neutral expectation as: E = Lxy — 
{L^ Ly), where L„ Ly, and are the relative lifespan of the Ax 
and Ay single knockouts and the Arc double knockout, respec- 
tively. Epistasis error was estimated as ffc = ^ff^^^ + aj^^ + a^^^ . 

Negative and positive interactions were defined at the 95% 
confidence interval. 



PLOS Genetics | www.plosgenetics.org 



13 



February 2014 | Volume 10 | Issue 2 | e1004168 



Systematic Genetic Analysis of Yeast Longevity 



The use of an additive neutral expectation results in a bias to 
positive epistasis (Figure SIO), as has been suggested for epistasis 
derived from cell growth based phenotypes [59]. We therefore use 

the conservative multiplicative model in all our epistasis analyses. 

Genes sharing similar epistatic-interaction patterns were clus- 
tered hierarchically (complete linkage) based on the correlation of 
raw e values (Spearman's rank correlation was used to account for 
strong oudiers at both sides of the distribution of e). Significant 
interactions and hierarchical clusters were used to generate an 
epistasis network, which was visualized with Cytoscape 2.8 [80]. 

Standard chronological lifespan assay 

For validation experiments using the standard CLS method 
[49], overnight cultures in non-bulfered SC (starting from three 
isolated colonies) were diluted to 0. 1 OD units into 25 mL fresh 
non-buffered SC medium and incubated at 30°C with shaking at 
220 rpm in 125 mL flasks covered with aluminum foil. Viability 
was measured by plating aging cells onto YPD-agar plates and 
monitoring colony-forming units starting from day 4 after 
inoculation, which was considered to be the initial survival (100%). 

Fluorescence microscopy 

Yeast cells expressing the CUPlpr-GFP-ATG8 construct were 
grown to mid-log phase {ODg(i(i=\) in SC medium. Expression of 
the fusion protein was induced with 10 mM CUSO4 for 2 h, 
followed by FM4-64 labeling at a final concentration of 10 jJ-g/ml 
for 20 min. Labeled cells were incubated 90 min in YPD medium 
to allow for dye uptake and then starved for 7 hours in SD-N 
medium to induce autophagy. Images were taken using a 
fluorescence microscope (Leica DM600 B) with a digital camera 
(Leica DFC420 C). 

Western blotting 

The autophagic degradation of Pgkl was monitored by 
immunoblotting of the Pgkl-GFP protein fusion as described 
[54]. Cells were grown overnight at 30°C to stationary' phase 
{ODgoo>A) in SC medium. 20 ODgog units of cells were harvested 
and starved for nitrogen in SD-N medium. At each time point 
after starvation, 2 ODgoo units of cells were lysed with 1.85M 
NaOH, 7.5% P-mercaptoethanol). Protein extracts were separated 
on SDS-polyacrylamide gels, transferred to a PVDF membrane, 
and probed with anti-GFP (Roche Applied Science 1 1814460001) 
and secondary antibody (Goat Anti-Mouse IgG coupled to HRP, 
31430 Thermo Scientific); reactive band were detected by 
chemiluminescence (GE Healthcare Life Sciences, RPN2209). 

Supporting Information 

Figure SI Competition-based assays provide rephcable CLS 
data. (A) Scatter plot comparing relative lifespan data (Z.) from two 
independent measurements of 720 deletion strains; r is the 
Pearson's hnear correlation coefficient. (B) Scatter plot comparing 
measurements obtained from two biological-replicates. WT and 
235 deletion strains (different from those shown in panel a) were 
independentiy tagged with RFP or CFP and the mutant's relative 
lifespan was measured in dye-swap experiments. 
(PDF) 

Figure S2 Competition provides increased quantitative resolu- 
tion for CLS measurements. The absolute half-life of WT and six 
gene-deletion strains was measured in rephcate [n = 24) by 
monitoring the outgrowth of either individual (gray box plots) or 
of competing populations (this study, yellow box plots). Outgrowth 
of individual cultures was monitored by the method of Murakami, 



et al. (2008, J. Gerontol. 63:113) adapted to our automated cell-assay 
station. For each outgrowth culture at a given age, viabilit\' was 
defined as the ODgoo „„ reached at a fixed point in time (10 hrs) 
and normalized to age zero (five days after inoculation); the half 
life of each strain was defined from the linear regression of 
ln(ODgoo nm) as a function of age (days). For competition-based 
assays, outgrowth was monitored by relative fluorescence as 
described and an absolute half life was obtained from s (see 
Materials and Methods). Starting from five days after inoculation, 
stationary-phas(; cultures were monitored for 27 days (individual 
cultures) or 21 days (tximpetition cultures). The dashed line 
indicates the median half life of the WT (21.7 days). 
(PDF) 

Figure S3 Relative survival ranks are insensitive to the choice of 
reference strain. (A) Heat map illustrates the relative survival 
coeflicients, s, of 24 RFP-labeled deletion strains aged in 
competition with 16 CFP-labeled strains that span the ^ scale. In 
the heat map, query and reference strains are ordered from short- 
lived (bottom left) to long-lived (top right); extreme negative and 
positive s values are thus concentrated at the top left and bottom 
right. (B) Cluster analysis of the Spearman's correlation coefli- 
cients (color bar) of the s values obtained for each CFP-labeled 
reference strain (treatment) against all RFP-labeled strains 
(observations). 
(PDF) 

Figure S4 Culture aeration has minor effect on CLS estimates. 

(A) CLS phenotypes at low and full a(-ration are highly correlated. 
Thirty-two mutant strains were aged in competition to the WT 
and outgrowth was monitored by relative fluorescence. Low 
aeration (horizontal axis) corresponds to 700 ^1 cultures in deep- 
well plates, as described (see Materials and Methods), while high 
aeration (vertical axis) corresponds to 8 ml cultures aged in test 
tubes shaken at 250 rpm and outgrown in microliter plates; r is the 
Pearson's linear correlation coefficient. (B) Semi-automated CLS 
profiling at high aeration confirms CLS phenotypes at low aera- 
tion. Plot shows the cumulative distributions of CLS phenotypes at 
high aeration for mutant strains defined as neutral, short hved, or 
long lived under the automated low aeration conditions used in 
this study. Cultures were aged in deep-well plates shaken at 
900 rpm on a Titramax vibratory shaker) Median L of both short- 
or long-lived strains was significantiy different from that of neutral 
mutant strains (/)<0.05, Wilcoxon rank sum test). 
(PDF) 

Figure S5 Comparison to growth-rate phenotypes allows 
identification of CLS-specific genes. (A) Short-lived mutant strains 
are more likely to be affected in their growth phenotype. Plot 
compares the distribution of growth-rate phenotypes of n(;utral 
(black dots) and short-lived mutant strains; p is the j6-value of the 
Wilcoxon rank sum test. (B) Venn diagrams show the numbers of 
mutant strains with CLS or growth-defect phenotypes; overlap 
indicates mutants affected in both CLS (short or long lived) and 
growth rate. CLS phenotypes were from this study (Table SI); Gis 
the relative growth rate, expressed as the average value of data 
obtained from Breslow d al. (2008, Mat. Metliois 5(8):711) and 
Costanzo et al. (2010, Science 327(5964):425). A mild/ strong grovrth 
defect (top) was defined as G<0.95, while a strong growth defect 
(bottom) was G<0.9. 
(PDF) 

Figure S6 Overlap of CLS and growth phenotypes increases only 
to a small extent under non-fermentative growth conditions. Scatter 
plots comparing relative lifespan (Z) in aging medium to growth rate 
(G) in (A) fermentative (dextrose) or (B) non-fermentative (ethanol-l- 
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glycerol) yeast-extract peptone medium for a random subset of 
knockout strains (« = 257). Red circles indicate strains with 
growth defects (G<0.9). Labels indicate gene mutations with 
overlapping CLS and growth phenotypes only under non- 
fermentative conditions (YPEG). (C) Pie charts show the fraction 
of growth-defect phenotypes under fermentative (YPD, left) or 
non-fermentative (YPEG) conditions for all tested mutant strains 
with a short- or long-lived CLS phenotype. Strains were growth 
in 150 (tl of YPD or YPEG medium in 96-weU plates shaken at 
900 rpm. Absolute growth rates (gr) were obtained from the 
linear fit of ln[ODg()g „J as a function of time; growth relative to 
the WT is G-grJgTM- 
(PDF) 

Figure S7 Our experimental approach can be readily adapted to 
measure lifespan independendy of potential regrowth. Cumulative 
frequencies of relative lifespan measured in water for mutant 
strains classified as short-lived (cross signs, n = 262) or long-hved 
(plus signs, n = 243). Plot shows strains that were confirmed in 
water as short-lived (magenta circles) and long-lived strains (cyan 
circles), based on the distribution of 57 WT-control replicates 
(95% CI). Strains were diluted into 800 |J,1 of buffered SC medium 
in semi-deep well plates and transferred to water five days after 
inoculation. Every three days, cultures were washed twice in water 
by microplate centrifugation. Cultures were monitored for relative 
survival during 21 days and relative lifespan was calculated as 
described (see Materials and Methods). 
(PDF) 

Figure S8 Characterization of single and double knockouts for 
epistasis analysis. (A). Plot shows the average relative lifespan (L^) 
of single autophagy mutants in which each mutant's Lx value 
(black circles) was derived from two to up to six independent single 
mutant replicates; error bars are the standard error. Single- 
knockout strains in the horizontal axis are ranked from largest to 
smallest L^. (B) Scatter plot compares L data from two marker- 
swap repUcates available for each double-gene knockout; r is the 
Pearson's linear correlation coefficient. 
(PDF) 

Figure S9 Most single mutants used for epistasis analysis are 
unaffected in their exponential-growth rates. (A) Strains were 
grown in 150 |J,1 of YPD medium in 96-well plates shaken at 
900 rpm. Absolute growth rates (gi) were obtained from the linear 
fit of IniODfioo nm) as a function of time during exponential phase. 
(B) Histogram shows the growth rates for sdl single knockouts 
relative to the WT, calculated as G = grjgr^rr. 
(PDF) 

Figure SIO Comparing the spectrum of epistasis derived from 
different neutral expectations. Epistasis was defined from (A) a 

multiplicative neutral expectation, &mulUpUcative = Lxy-{Lx-Ly), or 
(B) an additive neutral expectation, &addiuve = L^y — — Ly + 1. In 
each case, the spectrum of epistasis is shown for interactions 
among non-core-autophagy (blue histograms) and core-autophagy 
genes (purple histogram outiine, n = 75). Pie charts show the 
fraction of significant negative (e<0) and positive (e>0) interaction 
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